Theory for spin and orbital orderings in high temperature phases in YVO3 
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Motivated by the recent neutron diffraction experiment on YVO3, we consider a microscopic model 
where each l/"^"*" ion is occupied by two 3d electrons of parallel spins with two fold degenerate orbital 
configurations. The mean field classical solutions of the spin-orbital superexchange model predicts 
an antiferro-orbital ordering at a higher temperature followed by a C-type antiferromagnetic spin 
ordering at a lower temperature. Our results are qualitatively consistent with the observed orbital 
phase transition at ~ 200K and the spin phase transition at ~ 114K in YVOs. 



I. INTRODUCTION 

The transition metal perovskite oxides exhibit many 
interesting physical phenomena. In some of these com- 
pounds, the orbital degrees of freedom play an impor- 
tant role in their magnetic properties due to the strong 
spin-orbital coupling Q, |^. Examples include the 
Mott-Hubbard type insulators YVO3 and LaVO:^ , which 
show very unusual magnetic properties. Although the 
early experiments on YVO^ and LaVO^ were reported 
back in the mid 1970's 0, @i there has been re- 
newed interest in the past decade on these materials 
lllillESlilElElIllEi- There are two magnetic 
phases in YVO3: C-type antiferromagnetic order (fer- 
romagnetic chains along the z-axis which stagger within 
the x-y plane) at temperature 114K > T > 77K, and 
G-type antiferromagnetic order (staggered in all three 
directions) at temperature T < 77K 4, A H- The 
magnetic order in LaVOs is always C-type. The micro- 
scopic mechanism leading to the difference between these 
two compounds is still under investigation, and it might 
be related to the fact that at room temperature the cubic 
crystal structure is significantly distorted in YVO3 but 
almost undistorted in LaVO^. It is generally believed 
that the relevant orbital degrees of freedom, the degen- 
erate or almost degenerate 3d — t2g states are crucial to 
the observed magnetic properties. 

There have also been interesting theoretical studies re- 
lated to these magnetic behaviors 13, [13, H EES El 
ElEEE^- In particular, KhahuUin et.al considered 
I a spin-orbital Hamiltonian starting with 3-fold degener- 
' ate t2g orbitals, and compared the free energies between 
the C-type and G-type spin states in YVO3 by including 
an explicit Jahn- Teller ener gy i n the model. 

Very recently, Blake et.al [l^ reported neutron diffrac- 
tion experiment in YVO3, which shows clear evidence 
that the orbital ordering has a sudden change from high 
temperature G-type to low temperature C-type at the 
77K magnetic phase transition, manifested by a change 
in the Jahn- Teller type of distortion. The data also show 
clear evidence for the orbital transition from high tem- 



perature disordered phase to the G-type ordered phase at 
~ 200K. This has motivated us to study the spin-orbital 
ordering in YVO3. 

In this paper, we consider a microscopic model for in- 
sulating YVO3 , where each V-ion has two electrons with 
parallel spins favored by the large Coulomb repulsion and 
the Hund's coupling. The distortion in the cubic crystal 
structure already present at room temperature splits the 
degeneracy of the three t2g orbitals so that the d^y or- 
bital is favored by the crystal field. In our model, we take 
it as always singly occupied , while the other electron oc- 
cupies either the d^z or the dyz state. This description is 
consistent with the neutron diffraction experiment 
We consider the superexchange interaction of the model 
and derive an effective Hamiltonian for YVO3. We then 
study the mean field classical solutions of the model, and 
examine the spin and orbital orderings. We find a G- 
type orbital ordering at a higher temperature followed 
by an additional C-type spin ordering at a lower temper- 
ature. Our result is consistent with the observed orbital 
phase transition at ^ 200 K, and spin phase transition 
at 114K in YVO3. In this scenario, the orbital order- 
ing at ~ 200K is of the electronic origin, and the lattice 
distortion at ^ 200K observed in the experiment is a con- 
sequence of the orbital ordering and the electron-lattice 
coupling. The superexchange interaction alone consid- 
ered in our model does not explain the phase transition 
at 77K, which may require other interactions such as the 
Jahn- Teller effect as proposed in previous articles [TRIl^ . 

This paper is organized as follows. In Section II, we 
examine a multi-band Hubbard model at electron den- 
sity 2 electron per site, and consider the limit of large 
Coulomb repulsion and the large Hund's coupling. We 
then derive an effective Hamiltonian based on the su- 
perexchange mechanism. In Section HI, we discuss the 
mean field classical solutions of the model, and examine 
the phase diagram for the orbital and spin orderings. A 
brief summary is given in Section IV. 
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II. MODEL 

In YVO3, the vanadium electron configuration is 3cP. 
The compound has a cubic crystal structure, and each 
V ion is surrounded by six oxygen ions. Due to the cu- 
bic crystal field, the five-fold degenerate 3d orbitals are 
split into a higher energy doublet of Cg orbitals and a 
lower energy triplet t2g orbitals. At low temperatures 
and for low energy physics, the relevant orbitals are the 
three-fold t2g orbitals: d^y, dyz,dzx- In the strong cou- 
pling limit, the on-site Coulomb repulsion between the 
two electrons in the 3c? states and the Hund's coupling 
are much larger than the intersite electron hopping am- 
plitudes, the system is a Mott insulator with each V-ion 
having two localized electrons of parallel spins in two out 
of three degenerate t2g orbitals. This scenario appears to 
be consistent with experiments. 

As indicated in the recent diffraction experiment |l5l |. 
the cubic crystal is distorted at room temperature. As 
a result, the V-0 bond distances are anisotropic. Here 
we consider the structure at room temperature, where 
the V-0 bond distance along c-axis (perpendicular axis) 
is the smallest (see figure 1). This crystal structure fur- 
ther splits the t2g states. The dxy orbital has a lower 
energy, and becomes always singly occupied. The other 
d-electron is either in dyz or dzx orbital. In the diffraction 
experiment the data also indicate a smaller differ- 
ence in V-0 bond lengths in the xy plane, which we shall 
neglect here for simplicity. 
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FIG. 1: Idealized crystal structure for YVO3 studied in this 
paper. V — O bond distance along the 2-direction is shorter. 
The shown orbital represent configurations d^ydyz or d^yd^z- 

The atomic Hamiltonian j2l| is then given by i?o = 
Hi , where the sum over i runs all the V-sites, and 



(1) 



In the above Hamiltonian, c\^n„{cima) creates (anni- 
hilates) an electron of orbital m and spin a at site i, 
nima = cl^aC^ma- Ai = A2 = 0, and A3 = A < 0, 
with m = 1,2,3 representing orbitals dyz, dzx, dxy re- 
spectively. Umm' is the on-site direct interaction, and J 
is the exchange interaction, or the Hund's coupling. For 
the t2g orbitals, Umm = U = Umm' + 2 J for m' ^ m. In 
the case U,J > A, this Hamiltonian leads to an atomic 
ground state with each V — 3d^ ion having a total spin 
S — 1 with two-fold degenerate orbital configurations 
{dxy,dxz) and {dxy, dyz). This last restriction in orbital 
configurations is valid for YVO^ with strong lattice dis- 
tortion but not for LaVOs where the cubic structure is 
almost undistorted at room temperature. 

We next introduce the intersite hopping Hamiltonian 
Ht, given by 



(ij) mm' .a 



(2) 



where the sum runs over all the nearest neighbor V-V 
pairs, and t^'^ „j/ is the electron hopping integral between 
two sites i and j from orbital m to orbital m'. Since 
the most important contribution to the hopping inte- 
grals is from the path via the 2p state of the 0-ion be- 
tween the two neighboring V-ions, the hopping integrals 
are diagonal in the present problem due to the cubic 
symmetry. Namely, we have t^m,m' = tll^m.m'- There- 
fore, there are only two independent hopping parame- 
ters, t^i = t22 = t±, and t22 — ^33 = t^i — ^33 — 

with the super-index indicating the direction of the two 
sites. In the limit t±, ty «U, J, A, the system is an in- 
sulator with spin 1 on each V-ion. However, the virtual 
hopping introduces an effective intersite coupling of spins 
and the occupied orbitals. The effective Hamiltonian for 
H = Hq + Ht can be derived by applying perturbation 
theory to second order in t± or ty . 

Let \4>^^) = \sl,Tf,Sj,Tj) be a ground state of Hq 
for two V-ions where — 1,0, —1 is the spin z- 
component, and t is a pseudospin-1/2 operator for the 
orbitals: — 1/2 if dyz is occupied, and — —1/2 if 
dxz is occupied. The matrix elements between the unper- 
turbed ground states of the two V-ions can be calculated 
within the second order perturbation theory, and it is 
given by. 



min' ,aa' 



{<t)kl\Heff\4)i^ 



E 



{c^^i\Ht\I){I\Ht\^,^) 
Eq — Ej 



(3) 



where the sum is over all the intermediate eigenstates |/) 
of Hq corresponding to the eigen energy Ei, and Eq is 
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the ground state energy of Hq. Two-electron states with 
total spin S = 1 are given in ref. 21. The electronic 
configuration of the intermediate state |/) is 3(P on one 
V-ion and Sd} on the other. In the Appendix, we list 
all the states for V — and the corresponding energy 
difference Ej — Eq. The effective Hamiltonian can be 
derived from these matrix elements. Defining for each 
site a spin-1 operator S and a pseudospin- 1/2 operator 
T that act on the and degrees of freedom , it can 
then be expressed as below, 

^-// = E [I''i^:,'r,)S, ■ S, +L''(t„t,)] (4) 

where = x,y,z gives the direction of the bond (ij). 
The first term corresponds to spin-dependent orbital cou- 
plings while the second corresponds to orbital couplings 
which are spin independent. The first term also shows 
that the effective spin-spin couplings depend on orbital 
configuration. Equivalently, by defining = K"^ + K'^ 
and L" = K'^ — K'i , H^f / can be written as 

5][/^;(x„x,)(S..S, + l) (5) 

+X':(t.,x,)(S,.S,-1)] 

so that 2K'^ and 2K'^ are interpreted as the intersite 
orbital couplings for parallel spins (s| = 1, s| = 1) and 
antiparallel spins (sf = 1, s| = —1) respectively. We 
choose the energy unit to be i|/C/, and denote 77 — J/U , 
r,^ = 1/(1 - 377), 771 = 1/(1 + 277), and Q = t±/^. 
can be expressed in terms of parameters Q and 77, and 
are given below. 

if^^^=r73(r..r,,-i), 

if; = 2Q2^3(T,.f,-i), 

K^^)=a(r,,r,,-i) + J(l + 7,.) 

+ ^(l + ?7i)(^..+T,.), 

i^l=g2{2a(T,,r,,-i) + i(l + 77i) 

-\{r}i-^){rtr^ +r-T+) 

-\i^-rl^irtr^+rrT-)}. (6) 

In the above equations, a — — i(l + 2773 — 8771), 
and lx,y — —1 and -f-l respectively. For a bond in 
the z direction, where the dxy orbital is inert due to 
zero hopping amplitude and the dxz and dyz hopping 
is isotropic, our model is similar to the original Kugel- 
Khomskii model|j2^] with two differences. The first is the 



replacement of spin 1/2 by spin 1. The second is the effect 
of the dxy occupation which changes the Hund's coupling 
contribution to the on-site energies. 

We first discuss the intersite pseudospin couplings be- 
tween two parallel spins. In this case, the psuedospin 
has a SU (2) symmetry along z -direction. Along x- or 
y-direction, however, the virtual hopping integral for or- 
bital 2 or orbital 1 vanishes, so there is no exchange 
term in the pseudospin, and K^'^'^ is of the Ising form. 
The pseudospin coupling between the two V-ions of an- 
tiparallel spins is quite different. There is a linear term 
{Tiz + Tjz) along X- or 7/-direction, which either favors 
dzx or dyz orbital occupation to gain energy via the vir- 
tual hopping process. The pseudospin coupling along z- 
direction includes both the exchange term {t^t~ -\- h.c.) 

and the pair flip term {t^t^ + h.c). In spite of an 
isotropic matrix in the z-direction, the orbital Hamil- 
tonian is not SU(2) symmetric because of the presence 
of Hund's coupling. In particular, the pair flip term is 
related to superexchange processes involving those inter- 
mediate states listed in Table II that are split in energy 
due to Hund's coupling. To illustrate this point further, 
we consider a pair of V-ions along z-direction with an- 
tiparallel spins and pseudospins are r^^ = Tjz = 1/2. 
The relevant intermediate states in the superexchange 
are the states listed in the second and the fifth rows 
in Table II. Because these states have different energies, 
there is non-zero amplitude for the pseudospins to flip 
to Tiz = Tjz = —1/2. The pseudospin pair flip process 
is actually quite common in orbital physics. For exam- 
ple, there are pair flip terms in the effective Hamiltonian 
for spin- 1/2 systems with orbital degeneracy derived by 
Castellani et al. 21]. 

In the limit J/U = 0, 773 = 771 = 1, and a = 0, and 

so we have = Q^, and i^i'^'^^ = -%^(ri^ -I- J -I- 1. 
The orbital coupling between the two V-ions of antipar- 
allel spins vanishes, and the orbital coupling between the 
two V-ions of parallel spins remains to be pseudospin 
SU{2) symmetric along 2-direction and pseudospin Ising 
symmetric along x- or y-directions. For J/U = 0, the 
lack of the global pseudospin SU{2) symmetry is due to 
the anisotropic hopping integrals in the system. 

Our effective Hamiltonian here is similar to the Hamil- 
tonian proposed previously by KhaliuUin et al. with 
the following differences. These authors considered a 
model with 3-fold t2g orbital degeneracy, while we con- 
sider a 2-fold orbital degeneracy with dxy orbital being 
always singly occupied. Below we shall compare the two 
Hamiltonians by considering the the spin- orbital cou- 
pling between two V-ions along z-direction, namely Hij , 
with j = i + z. This may be carried out by imposing the 
orbital dxy being always singly occupied in the Hamilto- 
nian of KhaliuUin et al. We find that the TizTjz terms 
are the same in the two theories. However, the Hamilto- 
nian of KhaliuUin et al. does not include the pseudospin 
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flip term {t^t'^ + h.c). As we illustrated above, the 
pseudospin flip term is non-zero. As we shall see in the 
next section, this pair flip term does not affect the mean 
field results, which depend only on the z-components of 
pseudopsin in the present case. It will be an interesting 
question to examine if the pseudospin flip term is impor- 
tant to the orbital fluctuations. 



III. MEAN FIELD THEORY AND THE PHASE 
DIAGRAM 

We start with the classical solutions of i?e//- The 
Hamiltonian has a global SU(2)symmetry in spin space, 
so that we can assume the spin ordering along the z- 
direction. The Hamiltonian is invariant under the simul- 
taneous transformation of global Z{2) (reversing orbitals 
at all sites) and a 90" rotation of the lattice about the 
2;-axis. In general, we should consider orbital ordering 
along an arbitrary orientation. However, for the present 
problem, the orbital z-component terms are always larger 
or equal to the x- or y-component terms in K!^ of Eq. (6) . 
Therefore, we can discuss the classical solutions by con- 
sidering the z-component of the orbital ordering only \2l\ . 
In other words, the classical solutions are the same as the 
Ising solutions in the present case. In Table I, we show 
the energies per site in various classical states. Note that 
rji < 1, and rjs > 1, where the equality holds if and 
only if J = 0. We consider below the case with non-zero 
Hund's couphng J > 0. In this case, the two states listed 
in Table I with C-type antiferro-orbital (CO) configura- 
tion have higher energies. Also, we can see that the G- 
type antiferromagnetic spin (GS) and G-type antiferro- 
orbital (GO) phase has a higher energy than the CS-GO 
phase. Therefore, the ground state is either ferromag- 
netic spin (FS) and GO, or CS-GO. In both FS-GO and 
CS-GO phases, the orbital is antiparallel in all three di- 
rections, favored by the combination of the symmetries 
in hopping integrals (due to the cubic crystal symme- 
try) and the Hund's coupling. The asymmetry between 
the spin configuration along the z-axis and in the x-y 
plane is a result of the splitting of the dxy orbital level 
from the other two t2g orbitals. As expected, the FS-GO 
phase is energetically more favored at a larger J where 
the Hund's coupling dominates, and the CS-GO phase 
is more favored at a smaller J. It may be helpful to un- 
derstand these two possible ground states by examining 
the following limiting cases more explicitly. In the limit 
of large Hund's coupling, 77 — > 1/3, the terms in the en- 
ergy expression in Table I proportional to rja dominate. 
Hence the FS-GO phase has the lowest energy. In the 
limit J O"*", ryi — ^ 1 — O"*", and 773 ^ 1 -|- 0+, so that 
the CS-GO phase has the lowest energy. 

It should be noted that although the dxy orbital is al- 
ways singly occupied, the virtual hopping of the dxy elec- 
tron means that our model is not identical to a model 



with one electron occupying degenerate dxz and dyz 



or- 



bitals. In that model, for G-type orbital ordering, the 
ground state with non-zero Hund's coupling would al- 
ways be ferromagnetic. This is also true if we compare 
our model to a model with one electron occupying triply 
degenerate t2g orbitals |U 0| . 

We now discuss the finite temperature phases and their 
transitions. We introduce three types of order parame- 
ters, namely the spin order parameter m.i = (Siz), the 
orbital order parameter = {riz), and the spin-orbital 
order parameter qi = {SizTiz). We shall consider the or- 
der parameters corresponding to the FS-GO and CS-GO 
phases, since other ordered states are not energetically 
favorable. In both the FS-GO and CS-GO phases, we 
divide the lattice into two sublattices A and B accord- 
ingly. For the FS-GO ordering, we consider rrii = m for 
all the sites i, and Vi = r and qi = q for i at sublattice A 
and ri — — r and qi = —q for i at sublattice B. For the 
CS-GO ordering, we consider TOj ~ m, ri ~ r, and qi = q 
for i at sublattice A, and rrii — —m, ri = — r, and qi = q 
for i at sublattice B. We use a mean field theory to ex- 
amine the thermodynamically stable phases described by 
these order parameters, and neglect both quantum and 
thermal fluctuations. The effective Hamiltonian H^f / is 
then approximated by. 



Hmf = '^[aSiz + bTiz + cSizTiz + d) (7) 

i 

In the above equations, coefficient a, 6, c, d are functions 
oi rj , Q, as well as of the mean fields m, r, and q. They 
are given by, with the subscript upper (-) and lower (-I-) 
signs corresponding to the CS-GO and FS-GO phases, 
respectively, 

a = A^mr"^ + D^m, 

b — A^rn^r — B^r, 

c = -Azp(2m^r2 + q), 

2d = -C^m^r^ + E^m^ + Br^ - A^q"^, (8) 

where B = B-, and 

A± = 4(771 -f-a)(Q2±i), 
B± = 4(773 -a)(Q2^1), 

C± = 4(773 +a)(g2±i), (9) 
D± = (l+77i)(Q2±3)_(^^+a)(Q2±i)^ 
E± = -(1 + 77.)(Q2±^) + (773 +a)(Q2± 1). 

The mean field Hamiltonian can be solved easily to 
obtain the thermal averages of Siz, Tiz, and SizTiz, from 
which we obtain the following self-consistent equations 
for the order parameters m, r and q, (with (3 the inverse 
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TABLE I: Mean field classical ground state energies for various phases. 



Phase 


Energy per site 


G-type spin and G-type orbital (GS-GO) 


-|^5 + 2Q" + (1 + Q2)r?3 +3ryi j 


C-type spin and G-type orbital (CS-GO) 


-U5 + {1 + Q^)ri3 +3r;i ) 

o \ '2 2 J 


G-type spin and C-type orbital (GS-CO) 


-| ^10 + 30" + 2»7| + 3(2 + j 


C-type spin and C-type orbital (CS-CO) 


-f(^5 + r,3+3r,.) 


Ferro spin and G-type orbital (FS-GO) 


-2(1 + Q2)^3 



temperature) 



2 sinh(/3a) 
'(l + 2cosh(/3a)) 



1 



(r,,) = r = --tanh(/36/2) 



{SizTi 



sinh(/3c/2) 



(l + 2cosh(/3c/2)) 



(10) 



The free energy per site in the mean field theory is given 
by 



/ = -iln( 4cosh(/36/2)[l + 2cosh(/3c/2)] 



X [1 + 2 cosh(/3a)] 



(11) 



We solve the self-consistent equations for different or- 
dered states at various temperatures. The phases studied 
are (i). paramagnetic spin (PS) and para-orbital (PO) 
state (PS-PO) with m = r = q = 0: (ii). CS-PO state 
with C-type spin ordering 771 ^ and r ~ q ~ 0; (iii). 
PS-GO state with r 7^ and 777 = g = 0; (iv). CS-GO 
state; and (v). FS-GO state. In the states (iv) and (v), 
m,r,q ^ 0. When more than one set of MP solutions 
exist at a given temperature, we compare their free en- 
ergies to determine the thermodynamically stable phase 
diagram. 

In Fig. 2 and Fig. 3, we plot the phase diagrams 
obtained from the mean field theory for Q ~ 1.3 and 
Q — 0.5, respectively. The phase diagram for Q = 1 is 
qualitatively the same as that for Q = 1.3. The ground 
state is found to be CS-GO for smaller J/U, and FS-GO 
for larger J/U, consistent with our previous discussions. 
In general, the spin and orbital ordering occur at differ- 
ent temperatures. This feature is in distinction from the 
model for V2O3 for which spin and orbital order at the 
same temperature |2^ |2^ . Within the mean field the- 
ory, the phase transition between CS-GO and FS-GO is 
first order, and all other transitions between the different 
phases in Figs. 2 and 3 are second order. (The lattice 
distortion associated with the orbital ordering, which is 



PS-PO 
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CS-GO 



J/U 



FIG. 2: Phase diagram for the anisotropic hopping parameter 
Q — t±/t\i^ = 1.3 in the parameter space of temperature and 
relative Hund's coupling J/U . Temperature is shown in units 
of t y / ?7. Solid line is the phase boundary between para-orbital 
(PO) and G-type antiferro-orbital (GO) states, the dashed 
line is the boundary between paramagnetic (PS) and C-type 
antiferromagnetic (CS) states, the dotted line is the phase 
boundary between PS and ferromagnetic (FS) states. 



not included in our model, may change the nature of the 
phase transition.) 

In passing, we note that the spin-orbital ordering de- 
scribed by the order parameter q is always zero unless 
both the spins and the orbitals are ordered in the present 
system. This result indicates that the spin-orbital or- 
dering parameter q introduced in our mean field theory 
in addition to the spin order parameter 777 and the or- 
bital order parameter r may not be as significant in the 
present problem in altering the qualitative physics as in 
the S'C/(4) model [s^- 

In what follows, we shall focus on the phases relevant 
to YVOz, and discuss the sequential phase transitions 
from PS-PO to CS-GO. As we can see from Fig. 2, as 
the temperature decreases from the disordered state PO- 
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FIG. 3: Phase diagram for Q = t±/t\\ = 0.5. For details, see 
Fig. 2 caption. 



PS, the system first undergoes a transition to the G-type 
antiferro-orbital ordered phase. Only at lower temper- 
atures does the spin become C-type antiferromagnetic 
ordered. For smaller Q, the phase transitions depend 
on the Hund's coupling as we can see from Fig. 3. At 
intermediate J/J7, the orbital transition temperature is 
higher than the spin's, while at smaller J/J7, the spin 
transition temperature is higher than the orbital's. Since 
for YVO^, Q > 1, and the estimated values for U and 
J are [/ - A^.beV and J - QMeV [13, E3, our theory 
suggests orbital ordering first at a higher temperature 
followed by a subsequent spin ordering at a lower tem- 
perature in YVOz- This is qualitatively consistent with 
the experimental findings for YVO3 above 77K. 

As recently reported by Blake et.al the neutron 

diffraction experiment shows that the orbital ordering in 
YVO3 takes place at Too = 200K, which is far above the 
antiferromagnetic ordering temperature Tcs = 116K. 
The orbital ordering is evidenced by the changes of the V- 
O bond lengths in the xy-plane. Our theory is consistent 
with these observations. Orbital ordering and lattice dis- 
tortion are often observed simultaneously in experiments. 
It is usually difficult to distinguish if the lattice distor- 
tion is due to the orbital ordering or vice versa. Our 
theory suggests a scenario, in which the orbital ordering 
is of electronic origin, and the lattice distortion observed 
above 77K is a consequence of the orbital ordering and 
the electron-phonon coupling. 

In YVO3, as temperature decreases further, there is 
another phase transition at a lower temperature Tqs — 
77K, below which the system is in the G-type antiferro- 
magnetic and C-type antiferro-orbital state. There have 
been proposals [171 ll9j to attribute this lower temper- 



ature phase transition to the Jahn- Teller energy which 
favors C-type orbital ordering. It is an interesting issue 
to further understand the nature of the low temperature 
phase transition. 

IV. SUMMARY 

In summary, we have studied the electronic structure of 
the insulating YVO3, and derived an effective Hamilto- 
nian based on the superexchange interaction. We started 
with the atomic limit where each V — 3(P has a spin-1 
and two- fold degenerate orbital configurations {dxydxz) 
and idxydyz). This consideration is consistent with the 
recent neutron diffraction experiment at T > 77K. We 
studied the classical solutions of the model within mean 
field theory, and found G-type antiferro-orbital ordering 
at a higher temperature followed by a second phase tran- 
sition where the spins become C-type ordered. Our the- 
ory explains the orbital and spin ordering of YVO3 at 
temperatures T > 77K. While our model does not ex- 
plain the lower temperature phase of G-type spin order- 
ing, which may require considerations in addition to the 
superexchange interaction, our theory provides a starting 
point for understanding the unusual magnetic properties 

of yyOg. 

After we completed the present calculations, we learned 
of a very recent inelastic neutron scattering experiment 
of Ulrich et al. [U]. They reported an energy gap in 
the spin wave spectrum of YVO3 in the C-type spin 
ordering phase, and they interpreted it as an evidence 
for the orbital Peierls state along the z-direction. The 
classical solutions we study here do not predict any spin 
or orbital Peierls transition. Quantum fluctuations or 
electron-lattice interactions may be responsible for this 
unusual state [2(ll3ll|. 
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V. APPENDIX 

In this Appendix, we present all the intermediate eigen- 
states |/) and their corresponding energy differences with 
the ground state, Ej ~ Eq, of the the four t2g electrons 
in two V-ions. These states are used in the second order 
perturbation theory to derive the effective Hamiltonian 
Heff in the text. The atomic ground states of the sys- 
tem are 6x6 -fold degenerate, with each ion occupied 
by two electrons of parallel spins in the orbital configu- 
rations (dxydxz) or [dxydyz). The atomic ground state 
energy is E^ ~ 2{U — 3J). The excited states |/) are 
(y -M^)-{V- 3d^) with 6-fold degeneracy inV- 
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ion. Here we consider the limiting case L'^, J >> A, and 
neglect the effect of A on the V — 3(fi states. Within this 
approximation, these excited states are split into three 
multiples due to the Hund's coupling: 4x6 states with 
energy {U — 3J + Eq), 10 x 6 states with energy ([/ + Eo), 
and 6x6 states with energy ([/ + 2 J + Eq). The spin 
and orbital configurations of V — Sd"^ are listed in Table 
II. 
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TABLE II: Intermediate eigenstates |/) and corresponding energy differences Ei — Eq. The degeneracy shown is respected with 
the orbital configurations. 



Eigenstate 



(c||C2|C3|)|0) 

;y=(cJ|C2|C3| + cJ|C2|C3| + c||C2|C3| ) |0) 

,t J J 



4i4i)|o) 



— (c' c' 



^2t 



J-fct c""" - c' c' c' )\0) 



1|^2T^3T'' 



^(^T^i^J. 4j.4T4i)|Q) 



v^('^iT'^ii'^3T + '^2T '^21^31)!^) 
;yg(c||C|j^C3j^ + C2|C2|C3|)|0) 



Spin and orbital configuration 





^2 




1 
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Ei — Eo 



U 



U 

(3-fold) 



U 



U 



U + 2J 
(3-fold) 



